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Abstraet: A novel quantum imaging technique has recently been demonstrated in an experiment, 
where the photon used for illuminating an object is not detected; the image is obtained by interfering 
two beams, none of which ever interacts with the object. Here we present a detailed theoretical 
analysis of the experiment. We show that the object information is present only in the interference 
term and not in the individual intensities of the interfering beams. We also theoretically establish 
that the magnification of the imaging system depends on two wavelengths: the average wavelength 
of the photon that illuminates the object and the average wavelength of the photon that is detected. 
Our analysis affirms that the imaging process is based on the principle that quantum interference 
occurs when interferometric path information is unavailable. 


I. INTRODUCTION 

According to Bohr, comprehending the nature of a quan¬ 
tum system requires “a combined use of the contrast¬ 
ing pictures” of a classical particle and a classical wave 
[1]. Bohr’s complementarity principle [2] implies that the 
complete particle behavior and the complete wave behav¬ 
ior of a quantum system or entity are mutually exclusive. 
In other words, if a quantum entity behaves completely 
like a particle (wave) under certain experimental condi¬ 
tions, it does not display its wave (particle) behavior un¬ 
der the same conditions. To avoid confusion, we do not 
refer to a quantum entity as either “particle” or “wave”; 
instead, we use the term “quanton” (see, for example, 

[ 3 ]). 

The wave-particle duality can be illustrated by a 
lowest-order [4] interference experiment (e.g., Young’s 
double-slit experiment, Mach-Zehnder interferometer, 
etc.), in which a single quanton (e.g., photon, electron, 
etc.) is sent into a two-way interferometer (see, for ex¬ 
ample, [5]). If the quanton behaves completely like a 
particle, no interference can be observed at the output 
of the interferometer. It turns out that in this case it is 
possible to determine with complete certainty via which 
path the quanton has traversed. On the other hand, when 
there is absolutely no information on the path traveled 
by the quanton, perfect interference occurs—a behavior 
that characterizes waves. The relationship between in¬ 
terference and path information (wave-particle duality) 
has drawn the attention of several researchers (see, for 
example, [6-10]). 

The imaging process [11] of our interest is related to 
the wave-particle duality of photons. Let us consider two 
spatially separated identical light sources, 1 and 2, each 
of which has the ability of producing two photons at a 
time. These two photons are, in general, not identical 
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with each other and we label them by a and b. Suppose 
now that we select the a-photons from the both sources 
and send them into a two-arm interferometer under the 
following conditions: 1) photons from a particular source 
can travel through only one of the arms; 2) the sources 
emit at the same rate but in such a way that there is 
never more than one photon present in the interferome¬ 
ter at a time. In this case, although the a-photons are 
identical with each other, one can partially or fully ex¬ 
tract the interferometric path information by interacting 
with a 6-photon that is not sent into the interferometer. 
In such a situation, it is, therefore, possible to control the 
interference of a photon sent into the interferometer by 
using a photon that is not sent into the interferometer. 
This phenomenon has been experimentally demonstrated 
and discussed in Refs. [12, 13] and is often referred to as 
“induced coherence without induced emission”. 

The essence of our imaging technique [11] lies in the 
fact that the effect of interaction with 6-photons is ob¬ 
served in the first-order interference fringe pattern pro¬ 
duced by the a-photons. As for sources, we use two 
identical nonlinear crystals which generate photons by 
spontaneous parametric down-conversion. In Section II, 
we briefly recapitulate some basic results relating to the 
theory of spontaneous parametric down-conversion. In 
Section III, we then present a detailed analysis of the 
imaging method. Finally, we summarize our results in 
Section IV. 


II. ELEMENTS OF THE THEORY OF 
SPONTANEOUS PARAMETRIC 
DOWN-CONVERSION 

We mostly follow the theory of the process of sponta¬ 
neous parametric down-conversion developed by Hong 
and Mandel [14]. In this process a nonlinear crystal 
converts a photon (pump) into two photons (signal and 
idler) each of which has energy lower than that of the 
pump-photon. The combined energy of the signal and the 
idler photons is equal to the energy of the pump-photon. 
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When the pump beam is highly coherent and the down- 
conversion does not bring any observable change in the 
pump intensity, one can represent the pump by a clas¬ 
sical electric field Ep(r, t). In this case, the interaction 
Hamiltonian associated with the process of parametric 
down-conversion can be expressed in the interaction pic¬ 
ture as (cf. [13, 14]) 


= + H.C., 

( 1 ) 


where ^ represents the nonlinear electric susceptibility 

tensor of the crystal, E^~^(r, t) and Ej~^(r, t) are the 
negative frequency parts of the quantized electric fields 
associated with the signal and idler, respectively, D is 
the volume of the crystal, H.c. implies Hermitian conju¬ 
gation, and there is summation over the repeated indices 
q which label three mutually orthogonal directions 
in space. 

The pump, the signal and the idler fields may oscillate 
at different optical frequencies. In general, the suscepti¬ 
bility of the crystal depends on these frequencies. The 
Hamiltonian in Eq. (1) is therefore often expressed by de¬ 
composing the optical fields into several modes (see, for 
example, [14]). The positive frequency part of a quan¬ 
tized electric field inside the crystal can be represented 
by the expression [15] 


E*^+^ (r, t) = ^ a(k, a) exp [i(k • r — ujt)] ek,a S(k, a), 

k,cr 

( 2 ) 


where a = 1,2, labels two directions of polarization, uj 
is the frequency, k is the wave vector, ek,cr represents 
two generally complex, mutually ort hogonal unit vectors 
such that Gk • k = 0, (a(k, cr) = i^/huJ/{2eQn‘^{'k^ 
eo is the electric permittivity of free space, n(k, a) is the 
refractive index of the anisotropic, nonlinear crystal, 
is the quantization volume, and a(k, a) is the photon 
annihilation operator for the mode labeled by (k, a). Let 
us also decompose the pump field inside the crystal into 
plane wave modes and express it in the form 


Ep{r,t) = y] l/p(kp,(Tp)exp [i(kp ■ r - wpi)] ekp.ap- 

kp ,crp 

( 3 ) 


(cf. [16]) 


= j d^r EE {Ximgi‘^P^‘^s,‘^i)Vpi'kp,ap) 

^ kp,ap ks',as'k/jO-/ 

(ekp,<Tp); (eks,as)m (eL,aJ,«*(ks,crs)V(k/,(T7) 

exp [i{Lus +u)i - LOp)t] exp [z(kp - ks - k/) • r] 

aUk/jCr/)! + H.c., (4) 


where the subscripts P, S and / refer to pump, signal and 
idler, respectively. The quantum state of light generated 
by down-conversion at the crystal is given by the well 
known formula 


IV’(f')) = exp 



dt Hin{t) 


|vac), 


( 5 ) 


where |vac) is the vacuum state and t' is the interaction 
time. By expanding the exponential, Eq. (5) can be 
expressed in the form 


\ip{t’)) = jvac) + 



dt Hin{t) 


jvac) +- 


(6) 


III. IMAGING 

Let US consider a situation in which two identical nonlin¬ 
ear crystals NLl and NL2 are pumped by optical beams 
Pi and P 2 , respectively, generated by the same laser 
source (Eig. 1). The idler beam, Ji, generated by NLl 



FIG. 1: Illustrating the principle of the experiment. A laser 
beam (blue) is split into two beams Pi and P 2 which pump the 
nonlinear crystals NLl and NL2. The crystals produce signal 
(green) and idler (red) beams. The idler beam p is aligned 
with the idler beam p. The signal beams Si and S 2 are 
superposed by the beam splitter BS and one of the outputs of 
BS is detected by a camera. Dl, D2, D3 are dichroic mirrors. 


is transmitted through NL2 and is aligned with the idler 
The Hamiltonian, given by Eq. (1), now takes the form beam, p, generated by the latter. The signal beams Pi 
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and S 2 from the crystals NLl and NL2, respectively, are 
superposed by a beam-splitter, BS. One of the outputs 
of the beam-splitter is detected by an EMCCD camera. 

Suppose that the pump fields at the two crystals are 
given by the complex electric field vectors Ep^ ( 1 * 1 , ^ 2 ) and 
Ep 2 ( 1 * 2 , ^ 2 ), expanded in the form given by Eq. (3). Erom 
Eqs. (2), (3), (4) and (6) it follows that the quantum 
state of light generated by each individual crystal is given 
by the formula (cf. [17]) 


m')) 


t'D 


= |vac) + — E E E 


,aj. 


Vpj (kp,, crp.) a* (ks,, asj )a* (k/., ct/.) (ekp. ,ap .) ^ 
(eL, ,a,. ) exp [i{us^ + coi. - cop. )t'/2] 

sinc[(ws^- -u)p^)t'/2] 

exp [i(kp, - kp, - k/J • roj 

3 

I sinc[(kp^^ - ks^- - k/J„;„/2]| 

n=l 

|ksp<Ts,)s, 1+•••, (7) 


where j = 1,2 labels the two crystals, Ikp^., crp^. = 

(ks,, ) |vac)s^,, |kj., = a\^ (k/^, <7/.) |vac)^,, 

the volume integration has been carried out assuming the 
crystal to be a rectangular parallelepiped [18] of sides 
/i, / 2 , h with its center located at the point and 
sinc[x] = sinx/x; the sine terms lead to the two well 
known phase matching conditions associated with the 
process of spontaneous parametric down-conversion. 


A. Alignment of Idler Beams 

If the beam Ii is perfectly aligned with the beam / 2 , for 
each mode present in the quantized field of Ii there 
exists an equally populated mode in the quantized field 
E^^^ of / 2 . The perfect alignment of the idler beams 
can, therefore, be analytically expressed by the following 
formula: 


a/ 2 (k/,cr/) = a 7 ,(k/,cr/) exp[i0/(k/, a/)], (8) 

where the mode (k/,5/) is generated at NLl and is 
aligned with the mode (k/, aj) that is generated at NL2; 
(/)/(k/,a/) is a phase factor which can be interpreted as 
the phase gained by the mode (k/, 5/) due to propagation 
from NLl to NL2. 

The Hamiltonian given by Eq. (4) and, consequently, 
the state in Eq. (7) are expressed in quite general forms. 
We now simplify them using certain assumptions which 


are appropriate to our experiment. Let us assume that 
the signal and the idler can be treated as beams with 
uniform linear polarization both inside and outside of the 
crystals. In such a situation, we can drop the summation 
over the polarization indices ap.^ asj, (Ti. and can also 
write 

a(ks',,crs',) ~ as'(ws',), a{ki^,ai.) » (9) 

One can now replace the annihilation operator a(k, a) 
with a(k), the number state |k, a) with |k) and the sus¬ 
ceptibility tensor with a scalar quantity Eur- 

thermore, the alignment condition given by Eq. (8) re¬ 
duces to the form 

a/ 2 (k/) =a/i(k/) exp[i(t)i(kj)]. (10) 

It is to be noted that the relationship between k/ and k/ 
depends on the optical system used for aligning the idler 
beams. 



FIG. 2: Schematics of the imaging experiment. Positive lenses 
LI and L2 (both of focal length //), placed in the path of /i, 
form a 4-f system that images NLl on NL2. An identical 4-f 
system, consisting of positive lenses L3 and L4, is placed in 
the path of Si. A thin object is placed on a plane which is the 
back focal plane of LI and the front focal plane of L2. One 
of the outputs of BS is focused by a positive lens LO of focal 
length /o into a camera. NL2 and the back focal plane, F, of 
L4 are located at the front focal plane of LO. The lenses are 
assumed to be thin and ideal. The superposed signal beam 
is passed through a narrow-band filter, F, before entering the 
camera. 

In the experiment the two idler beams are aligned by 
the use of a 4-f lens system that images a central plane 
of NLl onto a central plane of NL2 (Fig. 2). We assume 
the idler beam axis to be along the optical axis of the lens 
system. A thin object which is intended to be imaged is 
placed at the back focal plane of the first positive lens, LI, 
of the 4-f system; this plane is also the front focal plane 
of the second positive lens, L2, of the same 4-f system. 
Clearly, the object is illuminated only by the idler beam 
that is generated by the first crystal. 

If a plane wave characterized by the wave vector k/ 
is incident on LI, it gets converted into a spherical wave 
that converges to a point , say, on the back focal plane 
of LI [Fig. 3(a)]. It then reemerges from this point as a 
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FIG. 3: Illustrating notations: (a) is a two-dimensional 
position vector lying on the object plane. The origin is chosen 
at the point O where the optical axis z meets the object plane. 
We neglect the limits due to diffraction, (b) The wave vectors 
kj and k/ are mirror images of each other with respect to the 
optical axis z. is the absolute value of the angle they make 
with the optical axis. 


diverging spherical wave. The amplitude of the diverging 
spherical wave can be determined from the amplitude of 
the incident wave and the complex transmission coeffi¬ 
cient object at point . The diverging 

spherical wave gets reconverted into a plane wave by the 
positive lens L2. This plane wave is characterized by a 
wave vector k/ which is different from k/, unless k/ is 
along the optical axis 2 : of the lens system. If one ne¬ 
glects the limits due to diffraction, one can say that a 
plane wave emerging from L2 contains information about 
one specific point of the object. 

Although the discussion of the previous paragraph is 
applicable to a classical field, it provides a guideline for 
treating the problem quantum mechanically. Since a 
quantized field is represented by decomposing it into sev¬ 
eral plane wave modes, one can say that one point on the 
object can transmit and refiect only one specific mode of 
the quantized idler field. Hence a single point on the ob¬ 
ject acts as a beam splitter only on one particular idler 
mode. Using the quantum mechanical treatment of a 
beam splitter ([19], sec. 12.12), one can now write the 
alignment condition [Eq. (10)] in the following form: 


S/2(k/) = 


{Pkj )«/i (k/) + ^' )ao (k/) 


X exp[i(/)/(k/)], 


( 11 ) 


where 3^ (p^^) is the transmission coefficient of the object 
at the point pk/, is the reflection coefficient at 

the same point when illuminated from the opposite direc¬ 
tion, So represents the vacuum field at the unused port 
of the beam splitter (a point on the object), ^/(k/) is the 
phase term mentioned earlier and = 1. It is 

evident that in absence of the object, i.e., when ^ = 1 
and = 0, Eq. (11) reduces to Eq. (10). 

If one neglects the limits due to diffraction, there is an 
one-to-one correspondence between k/ and k/, i.e., for 
every choice of k/ there is one and only one k/. When 
the focal lengths of LI and L2 are equal, k/ and k/ are 
mirror images of each other with respect to the optical 
axis 2 : [see Eig. 3(b)]. 


B. Detection System 


Let us now consider the detection system used in the ex¬ 
perimental setup. The signal beams generated by the two 
crystals are superposed by a 50 : 50 beam splitter (BS) 
and one of the outputs of the beam splitter is focused 
on an EMCCD camera by a positive lens LO (Eig. 2). 
A filter, E, is placed in front of the camera so that the 
light entering the camera has a narrow frequency band 
of mean frequency 



FIG. 4: Illustrating notations relating to the detection sys¬ 
tem: pk^ is a two-dimensional position vector lying on the 
image plane (camera screen). The origin is chosen at the 
point O' where the optical axis (also the beam axis) z' meets 
the image plane. The wave vector ks makes an angle Os with 
the optical axis. 


In absence of limits due to diffraction, LO would focus 
a plane wave with wave vector k^ at a point pkg, say, on 
the camera screen (Fig. 4). Hence the positive frequency 
part of the quantized field at the camera can be expressed 
as [20] 


E^s\p^s^t) 


(X exp 


-ius 





asi(ks) 


+ i exp 


iks ■ ro - iuJs 



r2(ks) ^ 


as2{ks), 

( 12 ) 


where Lj{ks) is the optical path traveled by the mode 
labeled by k 5 from the nonlinear crystal j to the cam¬ 
era (propagation inside the crystals has been neglected), 
|ks| = LOs/c and we have chosen ro^ = 0, roa — = rg. 


C. Formation of an Image 

For the sake of simplicity, we assume that the pump 
beams are well collimated, uniformly polarized and 
narrow-band with mean frequency up. In this case, 
the pump field at the crystal j can be represented by 
Vp. exp['i(kp • r — Opt)]. We choose the direction of kp 
to be along the direction of the optical axis. Using Eqs. 
(7), (9) and (11) one can show (see Appendix 1) that the 
quantum state of the field in this system can be approx- 
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imated by 

I vac) + ^ XIXI [9i^s^,i^h)Vp^ 

ks^ k/^ 

3 

{ ]^ sinc[(kp - ks^ - k/J„^„/2]} 

n=l 

sinc[(wsi +wp -u)p)t'/ 2 ] |ksi)s^ I 0 )g^ |0)o 

t'D r 

|^5(ws 2, W/JVP2 exp [i(kp - kp^ - k/J • ro] 

ks'2 kj2 

exp[-i()>/(k/2)]sinc[(wp2 +^^12 - i^p)t'/ 2 ] 

3 

{ sinc[(kp - ks2 - k/J„;„/2]} 

n=l 

+ '^'*(Pk.,)| 0 )s, | 0 )/Jk 5 .)s.|k/.)J], ( 13 ) 

where g{ujspt^i^ =X , w/,) (ws,) a|( w/,) exp [i {ujg^ 
uji- —(dp)t'/ 2 ], k/2 is the mirror image of k/2 with respect 
to the optical axis (beam axis) of the 4 -f system and 
we have suppressed the normalization coefficients. In 
the experiment, we choose crystals whose sides are ap¬ 
proximately 10“^m of length; since the wave vectors are 
characterized by corresponding optical wavelengths, the 
terms sinc[(kp^. — k^^. — k/^.)^/^/ 2 ] in Eq. ( 7 ) contribute 
only when kp^. — k^^. — kq. ^ 0 . This implies that 
the spatial phase matching condition kp^. k^^. + kq. 
holds to a very good accuracy. Furthermore, the term 
sinc[(co’5'^. + (jOj. — Co p)t'/ 2 ] leads to the temporal phase 
matching condition Cop ^ 00Sj + 001 .. 

The photon counting rate [ 21 ] at a point in the 
camera is given by 

7 ^(pk.) oc (^1 E^-\p^„t)E^+\p^„t) 1 ^), ( 14 ) 

where |T) is given by Eq. ( 13 ) and the quantized field is 
given by Eq. (12). It follows from a long but straightfor¬ 
ward calculation that apart from a proportionality con¬ 
stant 7 ^(pks) is given by 

l^(pkj « |Vpj2 + |Vpj2 

+ 2|VpJVp,||^(Pi^^)| 

cos 0s2(ks) - (/)Si(ks) - (t)i{ki) 

- arg[=^ {p^^ )] + ())p + kp • ro + Co , 

( 15 ) 

where k/ denotes a wave vector that is the mirror im¬ 
age of the wave vector k/ = kp — k^ with respect to 
the optical axis [Fig. 3 (b)], |kp| = Copjc, jk^l = Cos I 
|k/| = Coj/c, is the point on the object that is illu¬ 
minated by the idler mode k/, 05'^.(k^) = CosLj{'ks)/c, 


(l)p = arg[Vp 2 ] — arg[VpJ, and Cq includes all other phase 
terms. Equation (15) is the key equation of the theory 
of imaging. 

Let us first consider the situation in which no ob¬ 
ject is placed in the idler’s path, i.e., when |^| = 1 
and arg[^] = 0. We have already mentioned in Sec¬ 
tion IIIB that a point on the image plane (camera) has 
contribution only from one signal mode k^. Since the 
diameter of the signal beam cross-section in the cam¬ 
era is much smaller in dimension than the optical paths 
^^(k^) and the distance |ro| between the two crystals, 
the terms k^ • ro and ^^^(k^) can be treated as a slowly 
varying function of k^. Similarly, since the diameter 
of the first idler beam inside NL2 is much smaller than 
the distance between the two crystals, one can also ne¬ 
glect the k/ dependence of ^/(k/). These allow us to 
write 052 (^ 5 ') - ~ 0/(k/) ~ 0/o and 

k^ • ro + Co ~ Cq, where A^o, 4>io and Cq are constants. 
Equation (15) now reduces to 

^(PkJ«|Vpj 2 + |Vpj 2 

+ 2 |VpJ|Vp 2 | cos A^o — 070 + ffi Cq . 

(16) 

Since the right-hand side of Eq. (16) does not have any 
k^ dependence, it is clear that an almost uniformly il¬ 
luminated beam cross-section would be observed in the 
camera. The phase term 0p can be controlled in the ex¬ 
periment and by doing so one can modulate the intensity 
of the beam spot. It is evident from Eq. (16) that by 
changing the value of 0 p, one can achieve conditions both 
of constructive and of destructive interference, which are 
given by the following equations, respectively: 

As-o - 070 + 0PC + Cq = 2 A7r, (17a) 

As-o — 070 + 0PP + Cq = {2N + l)7r, (17b) 

where 0 pc' and <pPD are values of 0 p for constructive 
and destructive interferences [22], respectively, and N is 
an integer. Thus we have established the relations 

052(ks-) - (t^SiO^s) - 07(k/) + (ppc + Cq ~ 2N'k, 

(18a) 

052(k^) - 05i(ks') - 07(k/) + (jopD + Cq ~ {2N + l)7r. 

(18b) 

When the object is inserted in the idler’s path, it fol¬ 
lows from Eqs. (15) and (18) that the photon count¬ 
ing rates at a point in the camera under the conditions 
of constructive and destructive interference are, respec¬ 
tively, given by the formulas 

« |Vpj2 + |ypj2 

+ 21 VpJ I VpJ I ^(Pk, ) I cos(arg)]), 

(19a) 

« |Vpj 2 + |Vpj 2 

- 21 VpJ I Vp 2 11 ^(P£^ ) I cos(arg[^(P£^)]). 

(19b) 
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Equations (19) imply that an image of an absorptive ob¬ 
ject (arg[^(pkj)] = 0) and as well as of a phase object 
(|^(pk/)| = 1 ) would appear in the camera for both con¬ 
structive and destructive interferences (see [11], Fig. 3a). 
It further follows from Eqs. (19) that apart from a pro¬ 
portionality constant 

~ l=^(Pk,)|cos(arg[=^(pj^^)]), 

( 20 ) 

implying the background effect due to presence of the 
terms | Vp^ \ ^ and | Vp^ \ ^ can be eliminated by subtracting 
the photon counting rate obtained with destructive inter¬ 
ference from that obtained with constructive interference 
(see [11], Fig. 3d). It also follows from Eqs. (19) that 

7^(+)(pkJ +7^(-)(pkJ « 2 (VpJ2 + \Vp^\^). (21) 

This means that summing up the photon counting rates 
obtained by constructive and destructive interferences re¬ 
moves the image (see [11], Fig. 3c). Clearly, even if 
one uses an absorptive object that completely blocks the 
beam /i, the summation of the photon counting rates 
does not change. The fact that the information of the 
object appears only in the interference term shows that 
the imaging process is purely quantum mechanical in na¬ 
ture; this point is discussed later in further details. 

It is clear from the preceding discussion that a point 

in the object plane is imaged at the point pk^ in the 
image plane, where p^^ is the point at which a classical 

plane wave with wave vector k/ would be focused by the 
lens LI [see Fig. 3(a)] and pk^ is the point where the 
plane wave characterized by would be focused by LO. 
Since k/ is the mirror image of the wave vector k/ with 
respect to the optical axis [see Fig. 3(b)] and k/ is related 
to k^ by the phase matching condition k/ ^ kp — k^, the 
image that appears on the camera is not inverted [23]. 


D. Image Magnification 

So far we have neglected the effect of refraction at the 
crystal surface. However, one needs to consider this effect 
in order to obtain a correct value of the magnification. 
The wave vectors kp and k/ used thus far represent plane 
waves outside the crystals. The mirror image of k/ with 
respect to the optical axis z is k/ [Fig. 3(b)]. Suppose 
that the plane waves with wave vectors kp and k/ are 
represented inside the crystal by k|^ and kj, respectively. 
The phase matching condition can now be expressed as 

k'^ + k'/«k';, (22) 

where kp is the wave vector of the pump field inside the 
crystal. If k^ and kj make angles Og and 0'!^ respectively, 
with kp which is along the optical axis, it follows from 
Eq. (22) that 


cjpnp I sin 0^ \ ^ ujiui \ sin O'! \. (23) 


Let us choose two points on the object which are repre¬ 
sented by two-dimensional position vectors and p'~, . 

Suppose that their images at the camera are represented 
by the two-dimensional position vectors pk^ Pk' ’ 
respectively. The magnification is defined by the well 
known formula 


M = 


\p^s p'k'^ I 


(24) 


where its positive sign implies that the image is not in¬ 
verted. As already mentioned, the origins in the object 
and the image planes are chosen at the points O and O' 
where the corresponding optical axes (beam axes) meet 
the respective planes [see Figs. 3(a) and 4]. It readily 
follows from the theory presented in Section IIIC that 
O' is the image of O. By choosing the points and 

p^/ at O and O' respectively, we reduce Eq. (24) to the 
simplified form 


M = 1^. (25) 

\Pij\ 

If the signal plane wave k^, which is focused by the lens 
LI at the point pk^, makes an angle Os with the optical 
axis z' (Fig. 4), one has in the small-angle approximation 
IPksI = I/o Ian 6 ^ 5 ' I ^ Zolas' I- Similarly, one can show that 
l^kjl ^ where Oj is the angle made by the wave 

vector k/ with the optical axis [Fig. 3(a)]. It now follows 
from Eq. (25) that 


M = 


/o|^.g| 
fi\0i\ ■ 


(26) 


Since k^ and k/ are related to k.'^ and kj , respectively, 
by refraction at the crystal surface, using Snell’s law one 
obtains 


sin 6>5 = no sin 6 > 5 ', nj sinO'j = nosinOj, (27) 


where we have used the fact that k/ and k/ make the 
same angle with the optical axis [Fig. 3(b)] and the re¬ 
fractive index of air (no) has practically the same value 
for signal and idler. From Eqs. (23) and (27), it im¬ 
mediately follows that Os I sin 0s\ ~ ^7^/1 sin 0i\. In the 
small-angle limit, we thus obtain 

~ ^i\0i\. (28) 

From Eqs. (26) and (28) one finds that 

= (29) 

Jl^S fl^I 

Clearly, magnification of the imaging system depends on 
the ratio of the average wavelength of the signal to that 
of the idler. The dependence of the magnification on 
both wavelengths is a remarkable feature of our imaging 
process. 
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IV. CONCLUSION 


We have theoretically analyzed a recently demonstrated 
[11] quantum imaging technique. Although the experi¬ 
mental setup resembles an ordinary two-arm interferom¬ 
eter, the principle behind the imaging is purely quantum 
mechanical. If one imagines a classical two-arm inter¬ 
ferometer in one of whose arms an absorptive object is 
placed, both the interference term and the intensity con¬ 
tribution from the arm containing the object would de¬ 
pend on the transmissivity of the object. In our experi¬ 
ment, on the other hand, the intensity contribution from 
any of the crystals does not depend on the transmissiv¬ 
ity of the object. It is evident from Eqs. (19) that the 
information of the object is present only in the interfer¬ 
ence term. This also shows that the interference of signal 
beams is not due to the effect of induced emission (see 
also [12, 24]). This interference can only be explained by 
indistinguishability of the signal photon paths and hence 
the imaging process is directly related to wave-particle 
duality of photons. 

In this context, let us also have a close look at Eq. (13), 
which provides us with an expression for the quantum 
state |T) that has been used for explaining the imaging 
process. Since the pump source used in the experiment is 
a narrow-band laser, this state is obtained (see Appendix 
1) from the tensor product of the individual states, \ 2 pi) 
and I' 02 ) 5 generated by the crystals under the alignment- 
condition imposed by Eq. (11). This tensor product 
together with the alignment-condition implies the effect 
of induced emission. However, when the higher order 
terms present in |0i) and 102 ) can be neglected, the state 
|T) becomes identical with a state obtained by linear su¬ 
perposition of |0i) and 102 ) under the same alignment- 
condition (see Appendix 1). Since such a superposition is 
only allowed when there is no effect of induced emission, 
it is clear that photons generated by spontaneous para¬ 
metric down conversion play the key role in our imaging 
process [25]. 

Einally, although it is obvious, we would like to point 
out that the principle of our imaging also works if an en¬ 
tirely different lens system (or no lens system) is used in 
the experiment. A different lens system would only lead 
to a different value of the image magnification. Equation 
(29) shows that this magnification is equal to the product 
of two ratios: the ratio of focal lengths {fo/fi) and the 
ratio of wavelengths (A^/A/). The most remarkable fea¬ 
ture of this result is the presence of two mean wavelengths 
in the formula of magnification. This is a consequence of 
(a) the fact that the object is illuminated by photons of 
one mean wavelength while the camera detects photons 
of the other wavelength and (b) the phase-matching con¬ 
dition [see Eq. (23)]. The use of a different lens system in 
the setup might lead to a different expression for the im¬ 
age magnification; however, a dependence on two average 
wavelengths would always be present. 
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Appendix 1 

In this appendix, we illustrate the procedure used for ob¬ 
taining Eq. (13). We present a single-mode analysis, be¬ 
cause the procedure does not change when all the modes 
present in the quantized fields are considered. The single¬ 
mode and scalar version of Eq. (7) can be represented in 
the form 

IV^i) = |vac) + Gj |1)^^. \l)j. + ..., (30) 

where j = 1,2, the time dependence is suppressed, the 
coefficients are collectively denoted by Gj and the dots 
represent higher order terms containing higher powers of 
Gj. 

The single-mode version of Eq. (11) is given by 

Cii2 — (31) 

Erom Eq. (31), one immediately obtains that 

|1)/. = 11)/. I0)o |0),, |l)o] (32) 

Now using Eqs. (30) and (32), one can write 

|0i) = |vac) + Gi |1)^^ |1)/^ |0)^2 |0)q + ... (33a) 

IV’ 2 ) = |vac) + [^* lO)^^ \l)g^ |0)o 

+ £?'* |0)_5,^ |0)^^ |l)_g^ |l)o] +- 

(33b) 

Since both crystals are pumped by beams generated by 
a laser source (not a single photon source), the quantum 
state of light in the system is given by |T) = |0i) | 02 ), 
i.e., by 

1^-) = |vac) + Gl |l)_g^ |1)^^ |0)_g, |0)^^ 

+ G2e-*^^ [^* |0)sJl),Jl)5jO)o 

+ ^'* |0)_5.^ |0)j^ |1)_5.^ |1 )q] 

higher order terms, (34) 

where the higher order terms contains higher second or 
higher powers of the coefficients G. Since the rate of 
down conversion is very small, the these higher order 
terms can be neglected. Hence the state |T) can be ap¬ 
proximated by (neglecting the normalization coefficients) 

|T) ^ |vac) + Gl |1)^^ |l)/i | 0 )s 2 l ^)/2 
-h G2e 10 ) 5 ^ |l)7jl)5'2 l^)o 

+ 10 ) 5 , 10 ) 7 , 11 ) 5211 ) 0 ]- (35) 



Equation (13) is the multi-mode version of Eq. (35). 

It is to be noted that the form of |^) given by Eq. (35) 
can also be obtained by linear superposition of the states 
given by Eqs. (33a) and (33a), if one neglects the higher 
order terms. It is thus clear that if the experimental 


conditions are such that the contribution of these higher 
order terms is much smaller than that of the first order 
terms to a measurement of our interest, the measurement 
would yield the same result as in a case when the state 
1^) must be obtained by superposing \ 2 pi) and |' 02 )- 
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